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ABSTRACT 

We examine the evolution of inhomogeneities (fragments) of supernova ejecta in com- 
pact supernova remnants by means of hydrodynamical modeling and simplified ana- 
lytical calculations. Under the influence of intense post-shock cooling the fragments 
become strongly compressed as they traverse the hot shocked region between the re- 
verse and outer shocks of the remnant. We find that the most likely outcome of the 
interaction of fragments with the reverse shock and the hot shocked region is their 
disruption resulting in generation of secondary fragments. Secondary fragments arriv- 
ing at the thin and dense outer shell of the remnant give rise to brief X-ray flashes. 
Under suitable conditions the primary fragments may traverse the hot shocked region 
without being completely destroyed, to eventually reach the outer shell as dense, elon- 
gated structures. Collisions of such fragments with the shell are likely to give rise to 
powerful X-ray flares. 

Key words: supernova remnants - hydrodynamics - instabilities - shock waves - 
galaxies: active - X-rays: galaxies 



1 INTRODUCTION 

There is growing observational evidence that the matter ex- 
pelled in supernova explosions is not uniformly distributed, 
but instead it consists of clumps or fragments with a large 
variety of sizes, shapes and densities. For instance, the pres- 
ence of fragments has been detected, or at least suggested, 
in SN 1987A (Arnett, Fryxell & Miiller 1989, Lucy et al. 
1989, Hanuschik et al. 1993), SN 1993J (Spyromilio 1994, 
Wang & Hu 1994), and recently in Vela (Aschenbach, Egger 
& Triimper 1995, Strom et al. 1995). Fast moving metal-rich 
knots with velocities exceeding 5000 kms~^ have been iden- 
tified in Cas A (Braun, Gull & Perley 1987, Anderson et al. 
1994) and evidence for fragmented ejecta is also found in Ty- 
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cho (Seward, Gorenstein & Tucker 1983), Puppis A (Winkler 
et al. 1988), Kepler (Bandiera & van den Bergh 1991) and in 
a number of extragalactic remnants (Lasker & Golinowski 
1991, Fesen & Matonick 1993, Chugai 1993). This prepon- 
derance of dumpiness suggests that ejecta fragmentation is 
a common process, and its origin has been extensively stud- 
ied with the help of multi-dimensional numerical simulations 
(Miiller, Fryxell & Arnett 1991, Fryxell 1994 and references 
therein). It has been found that in a Type II supernova there 
are two layers, at the H/He and He/C-|-0 interfaces, which 
within minutes from core collapse develop strong hydrody- 
namical instabilities resulting in the formation of relatively 
dense fragments reaching, and probably exceeding, a factor 
of 10 in density contrast with the homogeneous interfrag- 
ment part of the ejecta. 

The effects caused by fragments in a SN ejecta were first 
explored by McKee (1983) and Hamilton (1985) by means 
of analytical and simplified numerical techniques assuming 
spherical symmetry (see the review by Franco et al. 1991). 
They found that the presence of fragments can increase the 
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thermalization radius of the ejecta, and that high column 
density fragments can move ahead of the main shock wave, 
driving a precursor which heats the ambient medium. 2-D 
numerical models of clumpy SNRs evolving in pre-existing 
wind cavities were studied by Tenorio-Tagle et al. (1991). 
The collisions between fragments and a pre-existing wind- 
driven shell can puncture the shell and induce a rapid mix- 
ing of different gas components. In addition, the propaga- 
tion of fragments across the remnant can be partially re- 
sponsible for the observed filamentary structures and X-ray 
halos in these regions. More recently. Franco et al. (1993a) 
explored, with an analytical approach and 2-D numerical 
simulations, the effects of clumpy ejecta in the evolution of 
multi-supernova remnants. As in the case of smaller bubbles, 
the fragments are likely to excite and puncture the remnant 
shell over time-scales of the order of 10^ yr, causing both 
temporal and spatial variations in the X-ray emission and 
introducing distortions in the shell structure. 

On the other hand, two new kinds of peculiar type II 
supernovae (SNe) have been observed recently, the very lu- 
minous type II radio supernovae and the so-called 'Seyfert 
1 impostors' (see review by Terlevich 1994). Both types 
seem to be associated with regions of active star formation. 
Seyfert 1 impostors are very bright in the optical contin- 
uum with very strong, broad Ha emission with no P-Cygni 
profile. Both the high Ha luminosity and the strong radio 
emission are interpreted as resulting from the interaction 
between the expanding SN ejecta and a dense circumstellar 
medium (Chevalier 1982, Terlevich et al. 1992, 1995, Chugai 
& Danziger 1994). 

Models of the interaction of SN ejecta with a high- 
density homogeneous circumstellar medium (CSM) match 
the observed features (spectrum, optical light curve, X ray 
luminosity and emission line widths) of these peculiar SNRs. 
Remnants evolving in a dense CSM (n > 10^ cm~'^) reach 
their maximum luminosity (L > 10^ L©) at small radii 
(-R < 0.1 pc) soon after the SN explosion [t < 20 yr) while 
still expanding at velocities of more than 1000 kms"^(Shull 
1980; Wheeler, Mazurek & Sivaramakrishnan 1980; Draine 
& Woods 1991; Terlevich et al. 1992). Because of their smaU 
dimensions, these remnants are designated 'compact' super- 
nova remnants (cSNRs) . A key feature of cSNRs is that ra- 
diative cooling becomes important well before the thermal- 
ization of the ejecta is complete, resulting in the remnant 
bypassing the quasi-adiabatic Sedov track. Consequently, 
the shocked matter undergoes a rapid condensation behind 
both the leading and the reverse shocks. Two concentric, 
high-density, fast-moving thin shells are then formed. These 
shells, along with the freely expanding ejecta and a section 
of the still dynamically unperturbed interstellar gas, are ir- 
radiated and ionized by the photon field produced by the 
radiative shocks. The resulting photoionized gas produces 
broad emission lines with line ratios that mimic those ob- 
served in the broad line region (BLR) of Seyfert galaxies 
(Terlevich et al. 1992; Franco et al. 1993b,c, Plewa 1995). 

Recent observational work by CoUura et al. (1994) in- 
dicates that the cSNR in the core of the nearby starburst 
galaxy M 82 may be varying at X-ray energies in less than 
a day. Collura et al. actually suggested that the variable 
source detected with ROSAT is associated with a massive 
binary, but examining the radio-maps of M 82 we found 
that its position coincides with that of the dominant radio 



supernova in M 82 (Kronberg, Biermann & Schwab 1985). 
Furthermore, Boiler, Fink & Schaeidt (1994) have discov- 
ered a luminous Starburst galaxy with an X-ray luminos- 
ity of ~ 10^^ ergs~^ that has variability on time-scales of 
less than a day and varies by more than a factor of two in 
amplitude. Inhomogeneities or fragments in the ejecta, we 
believe, are an important ingredient in modeling rapid vari- 
ability in cSNRs and starbursts. In this paper we focus on 
the evolution of clumpy (fragmented) SN ejecta evolving in 
a dense circumstellar medium characterized by the develop- 
ment of dynamical instabilities in the thin, cool and dense 
gaseous sheets formed as a result of intense cooling behind 
the shocks. In particular, we explore the possibility that the 
presence of fragments in the ejecta and/or in the ambient 
medium could lead to a rapid X-ray variability in cSNRs. 

Section 2 presents a set of 2-D numerical simulations 
following the fragments as they interact with the reverse 
shock and travel through the region of hot shocked gas be- 
tween reverse and forward shock fronts. Section 3 estimates 
the expected X-ray variability by means of analytical ex- 
pressions and Monte Carlo simulations. Finally, section 4 
summarizes our conclusions. 



2 NUMERICAL SIMULATIONS 

2.1 Methods, assumptions and input physics 

The hydrodynamic equations are integrated with the help 
of the AMRA code which combines the PPM method of 
Colella and Woodward (1984) with the AMR (Adaptive 
Mesh Refinement) approach first introduced by Berger & 
Oliger (1984) and subsequently developed by Berger & 
Colella (1989). The details of the method together with the 
results of test calculations are described elsewhere (Plewa 
& Miiller 1996). Briefly, the AMRA code operates on a 
series of nested meshes simultaneously, automatically cre- 
ating higher-resolution meshes and moving them over the 
base grid to follow the small-scale features of the model. A 
two-dimensional version of the code is used, able to work 
in spherical, cylindrical or Cartesian coordinates. For the 
present simulation spherical coordinates (r, 6, (j)) are chosen, 
implying a rotational symmetry with respect to the 6^ = 
axis. The basic grid extends over radii 0.0 < r < 1.3 x 10^^ 
cm and polar angles 0° < 9 < 10°, containing 64 and 20 
uniformly distributed points in r and 9, respectively. 

The code is allowed to generate up to three refined 
grids (1-3), each one consisting of several meshes. In the 
^-coordinate each refined grid provides a twofold resolution 
improvement with respect to the next coarser one, while in 
the r-coordinate the resolution improvement factors are 2, 3 
and 4 for grids 1, 2 and 3, respectively. Thus, the resolution 
of the finest grid is 8.5 x 10^^ cm in r and 6.25 x 10~^ deg in 6, 
equivalent to that of a uniform grid of 1536x160 points ex- 
tending over the whole computational domain. A new mesh 
is created whenever the local relative density contrast be- 
tween neighbouring cells exceeds 0.1 or the local relative 
pressure contrast exceeds 0.01 (the refinement procedure is 
not applied to the region occupied by the freely expand- 
ing ejecta). With that prescription for mesh generation a 
speedup factor of ~ 3.5 is achieved, meaning that a simula- 
tion carried out on a single grid of 1536x160 points would 
require 350 per cent of the CPU time actually used. 
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A reflecting boundary condition is applied at r = cm, 
6 = deg, and 6 = 10 deg, while a free outflow from the 
grid is allowed for at the r = 1.3 x 10^^ boundary. The equa- 
tion of state is that of an ideal monoatomic gas with solar 
composition. Complete ionization is assumed, resulting in a 
mean molecular weight /i = 0.615. Radiative, optically thin 
cooling of the shocked gas is allowed for, with the cooling 
function taken from Plewa (1995). For the total energy of the 
supernova explosion we use the standard value of 10"'^ erg, 
all of it in the form of kinetic energy of the ejecta. The explo- 
sion occurs in a uniform circumstellar medium (CSM) with 
a density ncsM = 10^ cm~"^ and a temperature Tcsm = 10^ 
K. 

Based on the results of Arnett (1988), we assume that 
the radial density distribution of the ejecta is well repre- 
sented by a p(r) oc r~'^ power law (the outer layer of the 
ejecta with a very steep density gradient is neglected, as it 
contains only a small fraction of the total mass of the ex- 
ploding star). We further assume that the velocity of the 
ejected gas is proportional to r, reaching about ten thou- 
sand kms~^ at the outer edge of the ejecta. As it is very 
likely that a typical massive supernova progenitor blows off 
most of its hydrogen envelope prior to the explosion, for the 
total mass of the ejecta a value of 5 Mq is adopted. 

As indicated in Section 1, the ejecta becomes frag- 
mented soon after the explosion due to Rayleigh- Taylor in- 
stabilities acting at the interfaces between layers with differ- 
ent chemical compositions (Arnett, Fryxell & Miiller 1989, 
Hachisu et al. 1992). Counting from outside, the first such 
interface is associated with the outer edge of the helium core, 
whose mass may vary between 4 Mq and 8 Mq depending 
on the initial (main sequence) mass of the progenitor. Thus, 
with a helium-core mass of ~4 M0 , the mass of the external, 
relatively smooth part of the ejecta in our model is equal to 
^1 Mq. The properties of the fragments are poorly known 
at present. In a simplified scenario one may describe the in- 
homogeneous part of the ejecta as an ensemble of cloudlets 
with various sizes, shapes, and densities immersed into an tn- 
terfragment medium (IFM). This IFM, apart from the large- 
scale radial density gradient, may be regarded as uniform. 
In such a scenario, and in the absence of internal motions, 
the fragments are frozen into the general diverging fiow, i.e. 
they expand homologously together with the IFM until they 
encounter the reverse shock of the remnant. According to 
Miiller, Fryxell & Arnett (1991), fragments as massive as 

0. 01 Mq can be found in the ejecta. In the present simula- 
tions a less extreme case is explored, with the mass of the 
fragment ranging from 3 X lO""* Mq to 3 x 10"^ Mq, while 
the whole ensemble of fragments with a broad distribution 
of masses is considered in the analytical part of the paper 
(Section 3). 

2.2 Initial conditions 

At t = the outer radius of the region occupied by the 
freely expanding ejecta is set equal to 3.25 x 10^^ cm (25 per 
cent of the base grid radius). Thus, our time count begins 
at about 2.5 yr after the explosion. For the first 3 yr of the 
evolutionary time the remnant is evolved without fragments, 

1. e. with the smooth IFM only. By then, the outer shock has 
propagated up to 7.8 x lO^'^ cm, while the reverse shock is 
located at 5 x 10^^ cm. Both shocks are still in the adiabatic 



evolutionary phase, and between them a hot shocked region 
extends, divided into two parts by a contact discontinuity 
which separates the shocked ejecta from the shocked CSM. 
We start the evolution of the fragments at t = 3 yr. At this 
time the mass of the shocked ejecta amounts to ~ 2.5 Mq, 
which means that not only the smooth hydrogen envelope, 
but also a fraction of the inhomogeneous ejecta (consisting 
of IFM and fragments) have been overtaken by the reverse 
shock. Thus, the fragments are inserted by the time when 
they are already interacting with the reverse shock, and we 
can follow the evolution of the shocked fragments and their 
possible interaction with the outer dense shell. 

In our simulations only one cylindrical fragment is fol- 
lowed. The fragment is generated at t = 3 yr in the freely 
expanding ejecta, with its front part almost touching the re- 
verse shock, and its shape approximated by a slightly elon- 
gated cylinder whose base is perpendicular to the radial di- 
rection (i.e. to the direction of motion). We define the e as 
the density contrast between the fragment and the inter- 
fragment media. The generation procedure simply increases 
the local density of the ejecta within the cylinder by a con- 
stant factor e, while the velocity and specific internal energy 
of the ejecta remain unchanged. The radius Rf and length 
l{ of the cylinder are the same for all simulations, amount- 
ing to 1.5 x 10^^ cm and 4.5 x 10^^ cm, respectively. With 
these dimensions, the range of fragment masses defined in 
the preceding section corresponds to an e-range of 3 to 30. 
The velocity of the front part of the fragment amounts to 
~ 4500 kms~^ (roughly equal to the largest velocity of the 
unshocked ejecta). 



2.3 Results 

Four different simulations have been performed, and we refer 
to them as cases A-D. Case A corresponds to (no fragment), 
and case B, C and D to e = 3, 10, and 30, respectively. In 
all cases the evolution is followed up to t = 11 yr, long past 
the thin shell formation phase. 

In case A, which serves as a reference case (not shown), 
the remnant evolves with a smooth ejecta without frag- 
ments. The shell begins to form at t ~ 6 yr at a dis- 
tance of ^ 2.5 X 10^^ from the outer shock front. Within 
less than 1 yr all the gas contained between the original 
shell formation site and the outer shock cools down and col- 
lapses into a very thin layer of shocked CSM (see Franco 
et al. 1994). This process occurs very rapidly and has been 
termed catastrophic cooling (see Plewa 1995 and references 
therein). Given the improved spatial resolution allowed by 
the adaptive method, this is the first time that the initia- 
tion of catastrophic cooling, and its subsequent migration as 
a cooling wave toward the outer shock, has been observed in 
2-dimensional simulations of supernova remnants. The cool- 
ing wave also propagates towards the contact discontinuity, 
condensing the shocked CSM onto the thin shell and causing 
the distance between the shell and the contact discontinuity 
to decrease. The hottest gas, which was the one shocked at 
the earliest moments, remains hot for a long time and even 
at the end of the simulation the remnant still contains an 
appreciable amount of high temperature gas. Our model res- 
olution provides an adequate qualitative description of the 
process, but most values and details cannot be completely 
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resolved. Thus, the present results only provide a reasonable 
guideline of the expected evolution. 

The basic evolutionary features of this case are also fol- 
lowed by the other three cases because the presence of a 
fragment in the flow does not alter significantly the global 
evolution of the remnant. In particular, the shell formation 
stage always proceeds in the same way and with the same 
time-scales. The main difference among models is in the de- 
tails of the fragment evolution, and in the emission and local 
structures generated at the sites of interaction. 

Snapshots of cases B-D taken at t = 4 yr are shown 
in Fig. |l| The stratified region extending up to 5.4 x lO'^^ 
cm from the center of the remnant contains the freely ex- 
panding (unshocked) ejecta. The sharp density contrast at 
r ~ 5.4 X 10^^ cm marks the momentary location of the 
reverse shock front. The hot shocked region extends to the 
right of the reverse shock, where a sharp and smooth right- 
hand boundary at r ~ 8.6 x 10^'' cm marks the location of 
the outer shock. The contact discontinuity separating the 
shocked ejecta from the shocked CSM is visible as a light- 
grey patch situated at r ~ 6.8 x 10^^ cm, roughly in the 
middle of the hot shocked region. The smooth indentation 
in the middle of the reverse shock, visible in all three panels 
of Fig. |l[ is the result of a perturbation caused by the pres- 
ence of the fragment (the shock slows down at the edges of 
the fragment). After the fragment has been entirely engulfed 
by the reverse shock, a conical self- reflect ion of the reverse 
shock is initiated, leading to a complicated pattern of sec- 
ondary shocks which can be best seen in the first panel of 
Fig. ^ and in all three panels of Fig. |^ (see also the detailed 
description in Klein et al. 1994). Secondary shocks are also 
driven into the rear side of the fragment. As all shocks sweep- 
ing through the fragment are strongly radiative (see Section 
3.1), the re-expansion phase described by Klein et al. (1994) 
is never reached. Instead, the fragment is compressed to the 
limits imposed by the resolution of the grid creating thin 
shells behind the penetrating shocks. Thus, aside from the 
Kelvin-Helmholtz instabilities described by Klein et al. , the 
shocked gas is also subjected to the thin shell instability 
and the obvious outcome of the interaction is the breakup 
of the original fragment into smaller and denser secondary 
fragments. At t = 4 yr (Fig. |l|) the breakup has only begun. 
In all three panels the fragment can still be recognized as 
a single object, but the effects of ablation (see also Kimura 
& Tosa 1991) and secondary fragment stripping are already 
visible. The main body of the fragment clearly shows the 
expected (see Section 3.1) dependence of the mean velocity 
on the initial density contrast (compare the positions of the 
fragment in all panels; compare also the strength of the bow 
shock in front of the fragment). The secondary fragments 
are very dense and move like solid particles inside the hot 
region, with the velocities of the fiow at the places and mo- 
ments of their creation (in particular, secondary fragments 
both overtaking the front edge of the primary fragment and 
lagging behind it may be seen in the lower two panels of Fig. 

0)- 

The subsequent evolutionary stages of case D are shown 
in Fig. ^ In all three panels only the magnified hot shocked 
region is displayed, bordered by the reverse shock (with a 
clear conical refiection pattern) on the left-hand side, and 
by the outer shock on the right-hand side. At t = 5 yr 
the formation of secondary fragments is completed. Being 



much denser than the shocked IFM, the fragments move es- 
sentially balistically across the hot shocked region. At later 
evolutionary stages, the effects of ablation and further frag- 
mentation can be seen. At the same time, however, some 
of the secondary fragments seem to be merging into larger 
entities (e.g. see the triple fragment at the symmetry axis 
and at r « 9 X lO^'' cm in Fig. ||). 

As mentioned above, between t = 6 yr (middle panel) 
and t = 7 yr (lower panel), a significant fraction of the 
shocked CSM behind the outer shock undergoes catastrophic 
cooling and collapses into a very thin and dense shell. Im- 
mediately after the dense shell has formed, the first frag- 
ments begin to collide with it. The results of the collisions 
are shown in Fig. ^ where the radiative activity of the mod- 
eled part of the remnant is illustrated. Obviously, given the 
large velocities involved (specially in the case of shocks faster 
than ~ 2000 kms~^) most of the thermalized kinetic energy 
is emitted in the form of X-rays (Plewa 1995), so that the 
X-ray and total light curves are very similar. 

The light curve of case A (Fig. m) is qualitatively much 
the same as the light curves obtained by Plewa (1995) for 
plane-parallel radiative shocks. For the first five years the 
luminosity of the remnant increases steadily due to the in- 
creasing cooling of the shocked CSM. At t = 5 yr the first 
mass elements of the shocked CSM reach the minimum of 
the cooling curve and enter the catastrophic cooling phase, 
in which the cooling rate increases with decreasing temper- 
ature. As a result, an X-ray/UV flare is generated, followed 
by a deep luminosity minimum between t = 5 yr and t = 6 
yr. The second flare at i = 7 yr is originated when the cool- 
ing wave arrives at the outer shock, and strong secondary 
shocks are driven into the cool CSM, completing the forma- 
tion of a thin and dense shell (see Plewa 1995 for a detailed 
description of this process). The differences in luminosities 
between the reference case A and the ones with fragments 
are illustrated in Fig. The first, low-amplitude peak is 
associated with the catastrophic cooling regime behind the 
shocks driven into the fragment. The total energy emitted in 
this peak is mainly dependent on the fragment density, with 
the largest luminosity corresponding to the largest e. The 
bombardment of the shell with secondary fragments begins 
at t « 6.25 yr (corresponding to the second significant peak 
in Fig. ^). As the distance between the formation site of 
the secondary fragments and the shell is equal to ~ 2 x 10^^ 
cm (see Fig. ^), one may estimate the speed of the fastest 
fragments arriving at the shell at ~ 3200 kms~^, which is 
below the initial fragment velocity (some 4500 kms~^). In 
contrast with the first peak, the total energies emitted in the 
second and subsequent peaks are correlated with the kinetic 
energy of the fragment impacting the shell. 

2.4 Discussion 

Our simulations show that an isolated fragment interacting 
with the reverse shock (and the hot shocked region) of a 
cSNR evolves through several distinct phases. Because all 
shocks are strongly radiative, the shocked fragment is com- 
pressed into a cool, dense, and unstable structure. Subse- 
quently, ablation and breakup of the unstable structure re- 
sults in the formation of secondary fragments with very high 
densities. Later, these secondary fragments collide with the 
outer shell of the remnant, giving rise to brief X-ray flashes. 
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Figure 1. Density distributions in cases B— D (top to bottom) at t = 4 yr. In case B (e = 3) the fragment is entirely sliocked, while the 
in case D (e = 30) the shocks driven into the fragment have swept about 75 per cent of its mass. The figure was obtained by reflecting 
all meshes with respect to the 6 = axis, i.e. the base grid extends over half of the area shown. Small asymmetries in plot details 
visible far away from the symmetry axis are artifacts of the plotting program. 
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Figure 2. The evolution of case D. Density distributions in the outer part of the remnant are shown at t = 5, 6 and 7 yr. The figure 
was obtained by reflecting all meshes with respect to the 9 = axis, i.e. the base grid extends over half of the area shown. Small 
asymmetries in plot details visible far away from the symmetry axis are artifacts of the plotting program. 
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Figure 3. (a) Light curve of the remnant with homogeneous ejecta, where Lgrid is the total energy radiated from the fraction of the 
remnant subtended by the grid. Luminosities for the whole remnant should be 132 times larger than Lgrid- (b) Differences in X-ray 
luminosities between remnants with homogeneous and inhomogeneous ejecta. Thin solid: e = 3, dotted: e = 10, heavy solid: e = 30. 



This general picture is not very sensitive to the fragment 
sizes, shapes, or initial density contrasts. It is conceivable, 
that at later evolutionary stages (when the distance between 
the outer shell and the reverse shock begins to decrease 
due to radiative cooling of the shocked ejecta), large-e frag- 
ments might collide directly with the shell before breakup 
(see Tenorio-Tagle 1994 for preliminary results of the evolu- 
tion with e = 100). These collisions should result in powerful 
X-ray flares. 

Our results provide an adequate qualitative description 
of the fragment evolution and interaction, but the details 
of the strong cooling phase and the thickness of the outer 
shell are not completely resolved. In fact, the shell may be 
more than two orders of magnitude thinner than our resolu- 
tion limit (see Section 3.2). Similarly, both the parameters 
of secondary fragments (specially sizes and densities) and 
their trajectories should be considered as rough, qualitative 
approximations. In general, an increase in numerical reso- 
lution would result in denser and thinner shocked layers. If 
the number of secondary fragments were not changing with 
increasing resolution, the characteristic time-scales of X-ray 
flashes would become shorter, and their amplitudes higher. 
Unfortunately, on the basis of the present simulations alone 
we are not able to decide if this assumption is realistic, i.e. if 
strong cooling combined with hydrodynamical instabilities 
results in formation of a few secondary fragments rather 
than in a complete destruction of the primary fragment. 
Nonetheless, it is expected that some fragments can collide 
with the shell before their complete destruction. The princi- 
pal conclusion of our simulations is that compact supernova 
remnants with fragmented ejecta should have a bombard- 
ment phase, in which dense clumps of cool gas collide at 
high velocities with the dense and cool outer shell creating 
X-ray flares. 



3 X-RAY VARIABILITY OF CSNRS WITH 
FRAGMENTED EJECTA 

The shell bombardment phase can last several years during 
which the remnant may be highly variable in the X-rays. 
Rapid X-ray variability is a direct observable diagnostic of 
the presence of fragments in cSNRs (see Cid Fernandes & 
Terlevich 1994 for the effects of the shell distortion on the 
emission line profiles). It is therefore important to estimate 
the variability properties one might expect to find in such 
objects. Since many different fragments may be involved si- 
multaneously in interactions giving rise to X-ray flashes, a 
direct numerical approach employing high resolution simu- 
lations is beyond our present capabilities. In the following, a 
statistical approach towards the X-ray variability of cSNRs 
will be applied, based on a simplified analytical expressions 
for the fragment evolution and the resulting interaction with 
the thin outer shell of the remnant. 

3.1 The evolution of fragments in the hot shocked 
region of the remnant 

Like in the preceding section, we assume that the inter- 
fragment gas is homogeneous. The thermal pressure of the 
shocked IFM can be evaluated from the conditions of mass 
and momentum conservation across the reverse shock front. 
If the IFM entering the shock has a density pifm and a neg- 
ligible pressure, and if it is processed by the reverse shock 
at a velocity Mifm, then the post-shock pressure is 

-Phsr = GmPifmttffm (1) 

where ^ifm = 4/3 for an adiabatic shock or 1 for an isother- 
mal shock. A rough, order of magnitude estimate for Wifm 
can be obtained using the planar two-streams approxima- 
tion, which yields 
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Olifm 
' 1 + difi 



(2) 



where Vcj is the velocity of the ejecta at the location of the 
reverse shock, au^i = (ncsMCcsM/nifm^ifm)^''^, and nitm is 
the number density of the interfragment medium. For an 
adiabatic shock with ncsM ~ nitm, the reverse shock pro- 
cesses the IFM at a speed of about 2wej/3, and the shocked 
gas streams away from the front at a speed of Uifm/4. 

PhBT defined by equation is the pressure that the 
fragments meet after passing through the reverse shock. 
From the jump conditions applied to a fragment that tra- 
verses the hot shocked region one easily obtains the strength 
of the lateral shocks driven into the lateral parts of the frag- 
ment: 



2 

fm^ifm 



For a fragment e times denser than the IFM, the relationship 
between Uf,i and unm becomes 



UC,l 



1/2 



el/2 



(3) 



As Uf^i < Uitm, the shocks processing the lateral walls of 
the fragment are weaker than the reverse shock propagat- 
ing through the IFM. The correspondingly lower post-shock 
temperature and higher post-shock density inside the frag- 
ments imply that the fragments cool faster than the IFM. 
In particular, high-e fragments may cool so rapidly that 
their internal shocks are practically isothermal (^f — 1). 
At the same time, since Uifm is typically of the order of sev- 
eral thousand kms~^, the reverse shock remains adiabatic 

tefm=4/3). 

To illustrate matters, let us consider the case of a reverse 
shock processing the IFM at a speed of 7000 kms~^, and 
assume that the density in the ejecta just ahead of the shock 
is lO'^ cm~^. Neglecting expansion effects, the cooling time 
for gas with solar abundances which enters a strong shock 
at a velocity Vs is (Franco et al. 1993b) : 



'^cool 



0.2(vs/n7) yr 
0.l{vi/m) yr 



if vs > 1.6 
if Vs < 1.6 



(4) 



where vs = Vb/W^ kms~^ and ny is the pre-shock density 
in units of 10^ cm""^, and the two regimes for «s larger or 
smaller than 1600 kms^^ are defined by the 'turning point' 
in the cooling function, at T ~ 3 x 10^ K. Hence, while 
the IFM takes about 1.4 yr to cool, an e = 100 fragment is 
processed by shocks of only ~ 600 kms~^ (using = 1 and 
Cifm = 4/3), and it radiates its post-shock thermal energy 
in just a few hours, while an e = 30 fragment would have 
Wf,i ~ 10^ kms~^ and a cooling time of a few days. 

Similar considerations apply to the shock which pro- 
cesses the front part of the fragment. Because of its forward 
motion through the hot shocked region, the front part of 
the fragment is exposed to an extra pressure Prci ~ PhsrWrob 
where Wrei is the relative velocity of the shocked fragment 
with respect to the shocked IFM, and phsr can be approx- 
imated by the IFM's post shock density (4pifm for ^ifm = 
4/3). The reverse shock therefore processes the front part of 
fragment at a speed 



1 + 



1/2 



(5) 



Vrcl 



The relative velocity Urci can be written as 

Mjfm _ UfJ 
Cifm (.i 

Replacing this relation into (^ we find that Ufj is between 
20 and 70 per cent larger than Uf,i for e between 10 and 
100. The resulting cooling times would therefore be similar 
to those derived for the lateral shocks (i.e., of the order of 
days or less for e> 30), still much shorter than the time the 
fragment needs to traverse the hot shocked region (~ 2 yr in 
our numerical models, see Section 2.3). Let us note that both 
the lateral and the front shocks are acting simultaneously, 
which means that the front shock propagates along the lat- 
eral walls through the gas that has been already shocked 
by lateral shocks, i.e. its velocity at the lateral walls of the 
fragment must be lower than the above estimate. Therefore, 
the above estimate for Ufj should be regarded as an upper 
limit for the effective (mass-weighted) front-shock velocity. 

Finally, v*{^ and vf, the velocities of the shocked IFM 
and the shocked fragment relative to the centre of the rem- 
nant are: 



"ifm 
5itm 



1 



(^) 



1/2 



°ifm 



where, for simplicity, we have assumed that ufj = Uf,i = U{. 
For flifin = 1, Cifm = 4/3, and = 1, these expressions yield 
^itm — ^cjfi and Wf = "Vcjll ~ (3e)~^'''^]. The fraction of 
the fragment kinetic energy which is thermalized and radi- 
ated away upon the interaction of the fragment with the 
hot shocked medium is 1 — (i^f /wcj)^. This loss amounts to 
only 10 to 30 per cent for e between 100 and 10, so that 
an appreciable energy is still available for the fragment-shell 
collision. 

Although the difference in lateral and front shock ve- 
locities is not too large (see the preceding paragraphs) the 
fragment is squeezed predominantly in the direction of mo- 
tion, becoming more and more fiattened. At the same time it 
is subjected to several instabilities which, as we know from 
the numerical simulations, result in the generation of sec- 
ondary fragments. Since low-e fragments may be completely 
destroyed and mixed with the shocked IFM, in the following 
we shall focus on fragments with e > 10, whose secondary 
fragments are very likely to collide with the outer shell at 
velocities of several thousand km s~^ (high-e primary frag- 
ments may even reach the outer shell before the disruption 
process is completed; see Section 2.4). 

An isothermal shock propagating through the fragment 
with a Mach number A4{ compresses the gas by a factor of 
Alf . Assuming that the sound speed in the unshocked frag- 
ment is 10 kms"'^, for the example with itifm — 7000 kms~^ 
discussed above one gets U{ « 600 kms~^ for e = 100, or 
10^ kms~^ for e = 30, with corresponding compression fac- 
tors between 3600 and 10000. Thus, fragments with a pre- 



shock density of ?s 
tion of motion) of 



10 cm and a length (in the direc- 
a lO'^* cm may reach densities ~ 10 



12 



10^^ cm ^, while their lengths decrease to 



lO^^-lO^^ cm 



(Fig. @a) 
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3.2 Fragment-shell collisions 

The cooled and condensed (primary or secondary) fragments 
eventually reach the outer shell of the remnant and collide 
with it at high velocities to produce hard X-ray flares. A 
collision between a fragment and the shell generates two 
new shock waves: a forward shock moving into the shell and 
a reverse shock moving into the fragment. The properties 
of these shocks, and their associated radiative bursts, can 
be estimated using the two-streams approximation. Again, 
assuming planar shocks, the velocities of the forward and 
reverse shocks are given by 



(l + a) 
a 



where V{ is the velocity at which the fragment impacts the 
shell, a = {us^is/nf^isy^^ and nt and Us are fragment and 
shell density, respectively. The shock crossing times for the 
fragment and the shell are respectively 



and 



' xm 



— = -^(l + a) — 

Vis 4fs «f 



1 (l + a) If 



V[ 



where l{ and Is represent the thickness of the fragment and 
the shell, respectively. Both forward and reverse shock waves 
remain radiative until one of them exits the high-density re- 
gions, i.e., until either the reverse shock overruns the frag- 
ment or the forward shock reaches the shell outer boundary. 
To compute the flare time-scale we shall approximate the 
duration of the collision (i.e., the time during which both 
shocks remain active) by the minimum of the two crossing 
times: Tcoii = niin(Txing, T"^ing)- Note that Tcoii is not the du- 
ration of the flare. The flare time-scales can be estimated as- 
suming that the radiative burst starts immediately after the 
shock passage and that the shocked gas remains radiative 
for a cooling time. In this case we have tI^,.^ ~ Tcoii -I- t^ooI 
for the fragment, and similarly for the shell. Seen by an ex- 
ternal observer, flares may be further stretched in time if the 
light travel time rnght = 2R{sm9/c between the near and 
far sides of the flaring regions is comparable to the crossing 
and cooling times. Here 6 is the angle at which the collision 
occurs, measured with respect to the observer's line of sight 
to the centre of the remnant, and -Rf is the characteristic size 
of the fragment in the direction parallel to the shell. For sim- 
plicity, we assume that the fragment approaching the shell 
has the shape of a flat cylinder, and we shall regard Rf as 
the radius of that cylinder. In general, the observed duration 
of the flare is 



(f 

ll+'^cool + flight 



(6) 



and the duration of the total flare is the maximum among 
r|jj^g and T|^^g. The thermal energy generated during the 
fragment- shell collision is quickly radiated away, and the to- 
tal energy contained in fragment and shell flares is 



Figure 4. Schematic representation of the evolution of a dense 
fragment inside a cSNR. (a) Interaction with the reverse shock 
and the hot shocked region, leading to the formation of a thin 
and dense pancake (for clarity, the lateral shocks are not shown), 
(b) Collision with the outer shell. Dark areas indicate shocked 
regions. 



-Eflaro — -E'flarc + -E'flarc — -EkiiiT 



■ niin(r?, 1) 



(7) 



where 77 = T^ing/Txing = als/li- The maximum flare energy 
corresponds to the case a 00, in which the fragment is 
completely halted at the shell and its kinetic energy is all 
radiated away. The spectral range in which the flare radiates 
is determined by the post-shock temperature, which for a 
shock velocity Vs can be calculated from 



1.2 



( V 

V10»cms-iy 



keV 



(8) 



To clarify the analytical model, the evolution of a frag- 
ment as it interacts with the reverse shock and later collides 
with the shell is schematically illustrated in Fig. ^. In this 
particular example, the fragment is completely shocked dur- 
ing the interaction with the outer shell, but since not all the 
kinetic energy is thermalized in the collision, the shocked re- 
gions keep moving after the end of the collision, producing 
perturbations in the surface of the shell. Smaller and/or less 
dense fragments may be completely stopped if the shell is 
sufficiently dense, in which case almost all the kinetic energy 
of the motion with respect to the shell is radiated away. 

To further illustrate the above considerations, let us 
compute the flare energy and time scale for a shocked frag- 
ment with density and thickness similar to those of the 



shell: Ic 



10 cm and nc 



where the shell parameters h and ris are taken from Ter- 
levich et al. (1995). The kinetic energy in the shell frame 
is Eii^ = 10"*^(mf/10"'')(i;f/10**)^ erg, with mf measured 
in M0 and Vf in cm s~^ (here we switch to characteristic 
fragment masses comparable to the Earth mass in order 
to account for the presence of secondary fragments). Since 
a = r] = 1, only half of this energy is actually radiated away 
during the flare (eq. The radius of a fragment with this 
length and density is -Rf = 8 x 10^^(mf/10"'^)^/^ cm, so 
the maximum light crossing time for the flare, correspond- 
ing to a collision at 6^ = 90°, is about 5200(mf /lO"^)^/^ s. 
The crossing and cooling times are identical for the reverse 
and forward shocks: Txing ~ 1500/(?;f /lO*) s and TcooI ~ 
40(wf/10®) s. Earth-mass fragments impacting the shell at 
speeds between 1000 and 10000 kms~^ would therefore pro- 
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duce flares with energies between 5 x 10^^ and 5 x 10''^ erg, 
and time-scales between 570 and 1540 s for collisions ob- 
served face-on, or 5770 to 6740 s for collisions observed edge- 
on. Fragments ten times more massive would produce flares 
ten times more energetic, whose durations would be in the 
same range as above for face-on collisions, but larger by a 
factor of about 3 for edge-on collisions due to the larger 
radii. The temperature of the shocked regions would be in 
the 0.5-50 keV range for V{ between 1000 and 10000 kms"\ 
We note that no strong X ray absorption is expected, at least 
for those collisions for which the cooling time is larger than 
the shock crossing time, since in this case the radiation goes 
through a column of highly ionized gas (see also Fig. 6b in 
Plewa 1995). 

It must be stressed that this simplifled scenario for 
fragment-shell interaction (thin and perfectly flat shocked 
fragments colliding frontally with a thin and perfectly flat 
outer shell) overestimates the efficiency at which the con- 
version of collision energy into a high-amplitude X-ray flare 
proceeds. Thus, the results obtained above actually provide 
a lower limit for the characteristic time scale and an up- 
per limit for the characteristic amplitude of X-ray flares in 
cSNRs, i.e., our estimates are based on the most favorable 
conditions to produce short and powerful flares. 

3.3 X-ray light curves and power spectra of 
fragmented cSNRs 

Both our numerical simulations and analytical results in- 
dicate that cSNRs may show substantial X-ray variability 
originating from the interactions between the compressed 
fragments of the ejecta and the dense outer shell of the rem- 
nant. X-ray light curves of cSNRs should therefore provide 
a straightforward diagnostic of the presence of fragments 
in such objects. In this section we explore the variability 
properties associated with fragmented cSNRs by means of 
Monte Carlo simulations. These simulations are important 
not only for comparison with future X-ray observations of 
cSNRs, but also to test their possible connection with ac- 
tive galactic nuclei (AGN). In the starburst model for AGN 
(Terlevich et al. 1995 and references therein) cSNRs are the 
source of high energy photons and the key ingredient of the 
broad line region. While the starburst model proved capa- 
ble of reproducing the main spectroscopic and photometric 
properties of broad line regions fairly well, it has not so 
far been able to explain the rapid X-ray flickering observed 
in AGN. Our simulations will thus be useful to determine 
whether fragment-shell collisions in cSNRs provide an ac- 
ceptable model for the rapid X-ray variability of AGN. In 
fact, given the lack of X-ray light curves for cSNRs, we shall 
concentrate the discussion on this possible application of the 
fragmented cSNR model. 

The short term cSNR light curve will be the result of 
a random superposition of many collision flares, associated 
with a multitude of fragments with a variety of masses and 
sizes interacting with the shell at different relative velocities 
and position angles. To compute light curve simulations we 
need prescriptions for (1) the time profile of each fiare, (2) 
the rate of fragment-shell collisions, and (3) the distribution 
of fragment properties (sizes, shapes and velocities) at the 
moment of collision. 

The time-profile of a given fiare can be calculated using 



the expressions derived above. We take full account of the 
light travel time effects, since the light crossing time across 
the lateral extent of large fragments can easily exceed the 
cooling and shock crossing times for collisions not seen face- 
on. The most critical assumption in the computation of the 
flare time-profile is the hypothesis that the collision is planar 
(see below) . By assuming a planar collision we are implicitly 
limiting our analysis to those sections of the fragments which 
do collide frontally with the shell. 

A rough estimate of the number of collisions per unit 
time {u) can be obtained dividing the total number of frag- 
ments by the duration of the shell bombardment phase 
(ivar), i.e., the phase during which the cSNR exhibits rapid 
X-ray variability. If a fraction of the ejected mass Mcj is 
in the form of fragments with an average mass mf, the mean 
rate of collisions is 



Cf A4j . (Cf/O.l) (Mej/lOMo) 



= 0.03- 



.r/yr) (mf/lO-«M0) 



(9) 



The distribution of fragment properties at the moment 
of collision with the outer shell depends on both the initial 
distribution and the details of the evolution of the ensem- 
ble of fragments in the shocked region of the cSNR. Since 
very little is known about the initial distribution of sizes 
and shapes, we are forced to adopt an ad hoc parameteri- 
zation. In the statistical simulations presented below we as- 
sume that masses and densities of the fragments as well as 
their radii R{ and velocities with respect to the shell have in- 
dependent power-law distributions. The shell thickness and 
density are kept fixed at 10^^ cm and 10^^ cm~^ respectively. 

Two of the many computed light curves are presented 
in Fig. U along with their corresponding power spectral dis- 
tribution (PSDs). Both models have fragment masses be- 
tween 10"* and 10"^ M©, radu between 10^^ and 10^^ cm, 
velocities between 1000 and 10000 kms~^ and densities be- 
tween lO^'^ and 10^^ cm~"^. The distribution of these prop- 
erties is however different in Figs. ^ and b: the logarithmic 
slopes (a's) of the mass, radius and velocity distributions 
are Qm = 0, Qr = and Ov = — 1 in model A (Fig. ^), and 
Qm = —1, ftr ~ —2 and Ov = in model B (Fig. BId), while 
fragment densities are homogeneously distributed (an = 0) 
in both runs. As a result of these distribution functions, 
model A favors larger and more massive fragments than 
model B. Accordingly, the ffares in A are more powerful 



(with an average energy iJfla 



3 X 10 erg) and longer 



5 X lO"*^ erg. 



("flare ~ 3 X 10 s) than in model B (iJfla 
"flare « 7 X 10^ s). The collision rate was adjusted to yield an 
output mean luminosity of the order of 10*^ ergs~^. The re- 
sulting mass in fragments impacting the shell per unit time is 
0.6 and 3.5 M0yr~^ in models A and B respectively. Clearly, 
collision rates as large as in model B cannot be sustained for 
a long time. (Note that the smooth background luminosity 
due to the homogeneous part of the cSNR is not included in 
the light curves. Such a background would have the effect of 
increasing the luminosity and diluting the variations.) 

Figs. ^ and b are examples of runs which produce light 
curves whose luminosities, r.m.s. variability and power spec- 
tra are similar to those observed in AGN. Rapidly vari- 
able AGN like NGC 4051, MCG 6-30-15 and NGC 5506 
have 'red-noise' PSDs with logarithmic slopes between —2 
and —1 in the 10~^-10~'^ Hz region (Lawrence et al. 1987, 
McHardy & Czerny 1987, Pounds & Turner 1987, McHardy 
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Figure 5. Left: Monte Carlo simulations of the X-ray light curve of fragmented cSNRs. The numbers in the plots are logu (in s~^, 
bottom-left), logL (in ergs~^, top-right) and the net r.m.s. variability (bottom-right). All light curves are normalized to the mean 
luminosity and binned in 500 s bins. Right: Power Spectral Distributions of the light curves. The power spectra were computed binning 
the periodogram using the technique of Papadakis &; Lawrence (1993). Power-law fits to the 10~^-10~'' Hz power spectrum are shown; 
the corresponding logarithmic slopes are listed in the bottom of the plots. 



1989, Green, McHardy & Lehto 1993, Lawrence & Papadakis 
1993). While some of the simulations have PSDs with log- 
arithmic slopes in this range (e.g. Figs. ^ and b), many 
others do not. 

As in any shot-noise model, the PSD is simply the av- 
erage of the power spectra of the individual flares (Lehto 
1989, Cid Fernandas 1995). Exponential flares with a decay 
time T, for instance, produce a PSD which is flat at low 
frequencies bending down to a law at high frequencies, 
the transition taking place within a decade around / = 1/r 
(Lehto 1989). Square pulses with a duration r produce the 
same effect, but with additional broad harmonics in the PSD 
at multiples of 1/r. Mathematically, in order to obtain a 
logarithmic slope between —1 and —2 over more than two 
decades in frequency we need a broad distribution of flare 
energies and time-scales such that slow flares dominate the 
low frequency part of the PSD while rapid flares dominate 
the high frequencies. In our physical model, this spread of 
energies and time-scales translates into fragments with dif- 
ferent masses, sizes, densities and velocities colliding with 
the shell at different angles. Getting the right PSD slope 
is a matter of finding a distribution of fragment properties 
which results in an appropriate balance of fiare energies and 
time-scales. Since many physical quantities are involved in 
the definitions of -Eflarc and rfiaro there are many combina- 
tions of the p{m{,V{, R!,nt) distribution function parame- 
ters which yield the same PSD slope, so it is difficult to con- 
strain the parameter space which yields an AGN-like PSD. 
In any case, the simulations show that at least for some dis- 



tribution functions the PSDs can be made similar to that of 
AGN. 

Though the fight curve simulations show that the model 
is in principle capable of reproducing the observed X-ray 
variability of AGN, the physical conditions required to ob- 
tain the sharpest (rflare ~ 1000 s) and more powerful bursts 
(Sflarc ~ 10^' erg) seem rather extreme. In order to reach 
such energies the fragments must have masses of 10~^- 
10~^ Mq and velocities of several thousand kms~^. These 
conditions, together with the requirement that the cool- 
ing and crossing times do not exceed 1000 s imply that 
the fragments should have radii of order 10^'' cm or more. 
Another important constraint is that the shell is very thin 
{Is ^ 10^^ cm), otherwise the shock crossing times would 
become exceedingly large. The strongest constraint, how- 
ever, is that the collision with the shell is almost perfectly 
planar. Slight distortions of the fragment can easily lead to 
fiare time-scales much larger than those resulting from a pla- 
nar collision. Our hydrodynamic simulations show that the 
fragments are far from being smooth, but the resolution we 
have been able to achieve is not sufficient for a definitive 
conclusion. 

Highly variable AGNs, as NGC 4051, MCG 6-30-15, 
and NGC 5506, pose yet another constraint to the model. 
In the present scenario it is clear that variations as large as 
those observed in these objects cannot be sustained through- 
out the cSNR lifetime, otherwise the mass and energy in 
the form of fragments would exceed those of the homoge- 
neous part of the remnant. The simplest solution to this 
problem would be to postulate that these objects are in an 
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evolutionary phase where rapid X-ray variabihty is at its 
peak: it either has been weaker in the past or will be so 
in the future. This is however in conflict with the results of 
Green, McHardy & Lehto (1993), who found no evidence for 
changes in the variability properties over intervals smaller 
than 3 years. It is nevertheless interesting to note that ob- 
jects like NGC 4151 and NGC 5548 do not show variations 
as large as NGC 4051, MCG 6-30-15 or NGC 5506, despite 
their similar X-ray luminosities (Green, McHardy & Lehto). 

To conclude, fragment-shell collisions are capable of 
generating 'chaotic' X-ray variability in cSNRs. With the 
free parameters suitably adjusted. X-ray light curves simi- 
lar to those observed in AGN can be obtained. Our feeling, 
however, is that in the case of extreme variability as that ob- 
served in AGN like NGC 4051, MCG 6-30-15 and NGC 5506 
the present model is strained too much. On the other hand, 
it seems definitely worthwhile to launch a program of X-ray 
monitoring of cSNRs, not only to test the hydrodynamic ef- 
fects predicted here, but to provide a strong empirical test 
of the possible connection between cSNRs and AGN. 



4 CONCLUSIONS 

We have studied the evolution of inhomogeneities (frag- 
ments) in supernova ejecta as the remnant interacts with 
a dense medium. Under such conditions, the ejected frag- 
ments, particularly the denser ones, experience a complex 
evolution and end up colliding with the remnant outer shell. 
The calculations reported in this paper show that: 

• Fragments are prone to further fragmentation inside 
the shocked region of the remnant. It is conceivable, how- 
ever that at least some fragments (those larger and denser) 
in advanced evolutionary phases of the remnant (when the 
distance between the forward and reverse shocks begins to 
decrease) may approach the outer shell as well defined enti- 
ties. 

• The shocked fragments eventually collide with the outer 
shell at velocities reaching several thousands kms^^. As 
both the fragments and the shell are already strongly com- 
pressed due to intense cooling, their densities at the moment 
of collision are very high (up to some 10^^ cm~^), and the 
shocks generated by the collision are strongly radiative. The 
outcome of such an energetic collision is an X-ray flare. 

• The observed total energy of the flare is essentially de- 
termined by the kinetic energy of the fragment in the refer- 
ence frame of the shell. The observed duration of the flare 
depends both on shell and fragment geometries and on the 
inclination of the collision plane with respect to the ob- 
server's line of sight. Under the most favorable conditions 
(a head-on collision of an entirely flat fragment with an en- 
tirely flat shell, collision plane perpendicular to the line of 
sight) the coUision energy can be thermalized and radiated 
in a fraction of an hour. 

• Fragment-shell collisions occurring randomly in space 
and time are in principle capable of generating light curves 
and power spectra similar to those observed in AGNs, but 
the conditions required to reproduce the sharpest and most 
powerful bursts seem contrived. 

On the theoretical side, future research should explore 
modifications of the basic scenario outlined in the present 



paper. For example, it is conceivable that X-ray flares can 
also be produced by collisions between ejected fragments 
and circumstellar clouds (see McKee 1983). The external 
cloudlets entering the hot shocked region would evolve simi- 
larly to the ejecta fragments, and a cloud-fragment collision 
should be similar to a fragment-shell collision. Another area 
for exploration is the interaction between ejecta fragments 
and the reverse shock thin shell. In the vicinity of the reverse 
shock, given that even low-e fragments are not yet severely 
distorted, the probability of a cloud-fragment collision is in- 
creased. Also, because of the resulting larger relative veloc- 
ities, one would expect stronger radiation flashes and with 
a shorter duration than those from collisions of ejecta frag- 
ments with the main shell. 

Further simulations with higher resolution will be neces- 
sary to uncover the details of fragment-shell collision. Length 
scales of the order of lO'^^ and 10^* cm are needed to pro- 
duce sharp X-ray flares. The corrugation of the shell in Fig. 
^has a characteristic wavelength of ~ 3 x 10^® cm), which is 
the expected value for the fastest growing mode of the thin 
shell instability in our simulation (the fastest growth rate 
appears on scales comparable to the thickness of the shell; 
Vishniac 1993). Thus, one can expect that an increase in 
numerical resolution would lead to a shell which is thinner 
and denser (see Sect 2.4), and with a finer corrugation scale. 
With our present computational resources, the required res- 
olution cannot be achieved but, given the rapid improve- 
ment in speed and memory resources, the problem could be 
approached in a few years. 

On the observational side, the basic verifiable predic- 
tion of the fragmented ejecta model is that cSNRs by them- 
selves should exhibit stochastic and rapid X-ray variability. 
Unfortunately, this prediction is also hard to test because 
bona fide cSNRs are hard to come by, being generated at 
unpredictable times and places. To date, the only relevant 
observation is that of the cSNR at the core of M 82 (Col- 
lura et al. 1994), which indeed show variations in scales of a 
day. Ideally, one would like to monitor a few cSNRs for a few 
years, but this would require large amounts of satellite time. 
SN 1988Z (Stathakis & Sadler, Turatto et al. 1993) recently 
observed as a strong X-ray source by ROSAT (Fabian & 
Terlevich 1996) is an obvious first target for such a project. 
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